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ABSTRACT 

Models of magnetically driven accretion and outflows reproduce many observational prop¬ 
erties of T Tauri stars, but the picture is much less clear for the Herbig Ae/Be stars, due to 
the poor knowledge of their magnetic field strength and topology. The Herbig Ae star PDS 2 
was previously included in two magnetic studies based on low-resolution spectropolarimetric 
observations. Only in one of these studies the presence of a weak mean longitudinal magnetic 
field was reported. In the present study, for the first time, high-resolution HARPS spectropo¬ 
larimetric observations of PDS 2 are used to investigate the presence of a magnetic field. A 
firm detection of a weak longitudinal magnetic field is achieved using the multi-line singular 
value decomposition method for Stokes profile reconstruction ((Bz) = 33 + 5 G). To gain bet¬ 
ter knowledge of typical magnetic field strengths in late Herbig Be and Herbig Ae stars, we 
compiled previous magnetic field measurements, revealing that only very few stars have fields 
stronger than 200 G, and half of the sample possesses fields of about 100 G and less. These 
results challenge our current understanding of the magnetospheric accretion in intermediate- 
mass pre-main sequence stars as they indicate that the magnetic fields in Herbig Ae/Be stars 
are by far weaker than those measured in T Tauri stars. 

Key words: stars: pre-main sequence — stars: individual (PDS 2) — stars: magnetic field — 
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1 INTRODUCTION 

Recent observations of the disk properties of intermediate mass 
Herbig Ae and late Herbig Be stars suggest a close parallel to 
T Tauri stars, revealing the same size range of disks, similar opti¬ 
cal surface brightness and similar structure consisting of inner dark 
disk and a bright ring. Models of magnetically driven accretion and 
outflows successfully reproduce many observational properties of 
low-mass pre-main sequence stars, the classical T Tauri stars. How¬ 
ever, due to the very poor knowledge of the magnetic field strength 
and magnetic field topology in Herbig stars, current theories are 
not able to present a consistent scenario of how the magnetic fields 
in Herbig Ae/Be stars are generated and how these fields interact 
with the circumstellar environment, consisting of a combination of 
accretion disk, magnetosphere, and disk-wind region or jets. 

As of today, only about 20 late Herbig Be and Herbig Ae stars 
have been reported to possess large-scale organ ized magnetic fields 
(e.g., iHubrig et alj|2009l : lAleciMi et al.ll2013^ . using low- and/or 
high-resolution spectropolarimetric observations. M oreover, only 
for the two Herbig Ae stars HD 101412 and V380 Ori jHubrig et all 
I 2 OI lal : lAlecian et al.ll2009ll . the magnetic field geometry has been 
constrained in previous studies. It is very likely that the rather 
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low detection rate of magnetic fields in Herbig Ae stars can be 
explained by the weakness of these fields and/or by rather large 
measurement uncertainties. Indeed, in the currently largest high- 
resoluti on spectropolarimetric survey of the magnetic field in these 
stars bv lAlecian et al.l ( l2013ah with 132 measurements for 70 Her¬ 
big stars, the measurement uncertainty is worse than 200 G for 35% 
of the measurements, and for 32% of the measurements it is be¬ 
tween 100 and 200 G, i.e. only 33% of the measurements showed 
a measurement accuracy below 100 G. Therefore, any new detec¬ 
tion/confirmation of the presence of a magnetic field in a Her¬ 
big Ae/Be star is important to increase the sample of magnetic 
Herbig stars. Clearly, only after we gain knowledge of the mag¬ 
netic field strength and its topology, can we start to understand how 
it interacts with the circumstellar environment. 

Polarimetric observations of the Herbig Ae star PDS 2 
(CD-53° 251) were obtained on two different epochs using the 
HARPS polarimeter. This star has been ide ntified as a Herbig Ae 
candid ate star in the Pico dos Dias Survey bv iGregorio-Hetem et al.l 
l ll992h . PDS 2 displays Ha in emission, and the pre-main sequence 
nature is indic ated by its IRAS colour s. The mass accretion rate was 
determined by [Pogodin et alj ( 12012h on eight epochs using differ¬ 
ent spectral accretion indicators from the near-UV to the near-IR. 
This work revealed short-term night-to-night variability, as well as 
long-term variability on a time scale of tens of days. 
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Using the 0.6 m R.E.M. telescope on L a Silla, a preliminary 
analysis of PDS 2 by iBemabei et al.l l l2007h showed the presence 
of 6 Scuti-like pulsations with three frequencies on a time scale of 
1.0-1. 7 h. The refined analysis of the same data by iMarconi et al.l 
( l2010h revealed pulsations at four frequencies with pulsation pe¬ 
riods between 0.9 and 2h. The relatively low projected rotation 
velocity, v sin i = 12 ± 2kms“*, and the fundamental parame- 
ter s, Teff = 6500K a nd log g = 3.5, were recently determined 
bv ICowlev et alj l l2014ll using high-resolution HARPS spectra. Ac¬ 
cording to their results, PDS 2 belongs to the class of Herbig Ae 
stars resembling the chemical composition of A Boo stars. 

Due to the relative faintness of PDS 2 (V = 10.8), the presence 
of a magnetic field was not previously studied using high-resolution 
polarimetric spectra. In the past, this star was included in two mag¬ 
netic studies based on low-resolution spectropolarimetric observa¬ 
tions obtained with FORS 1 (FOcal Reducer low dispersion Spec¬ 
trograph) mounted at the 8-m Kueyen telescope of the VLT. While 
in the first study by I Wade et al.l l l2007h no convincing evidence for 
the pre sence of a magnetic field in PDS 2 was found, iHubrig et al.l 
ll2009h reported on the detection of a mean lo ngitudinal magnetic 
field (Sj) = 103 ± 29 G. lBagnulo et al.l | |2012|) disputed this detec¬ 
tion after re-examining the FORS 1 data using their semi-automatic 
measurement procedure. 

In Sect. |2l we describe the observations and data reduction, 
and in Sect. [3] we discuss the methods and results of our magnetic 
field measurements. The density distribution of longitudinal mag¬ 
netic fields in magnetic late Herbig Be and Herbig Ae stars is pre¬ 
sented in Sect, [d] Finally, in Sect. |5] we discuss the significance 
of the obtained results for the improvement of our knowledge of 
the role of magnetic fields in intermediate-mass pre-main sequence 
stars. 


2 OBSERVATIONS AND DATA REDUCTION 

Two spectropolarimetric observations of PDS 2 were o btained with 
the HARPS polarimeter IHARPSpoh lSnik et alj|200^ attached to 
FSO’s 3.6 m telescope (La Silla, Chile) on the nights of 2012 July 
15 and 18. These spectra were originally recorded as part of FSO 
programme 187.D-0917(C) (PI: Alecian) and were downloaded 
from the FSO archive under request MSCHOFLLFR 77411. Fach 
observation was split into eight subexposures with an exposure time 
of 15 min, obtained with different orientations of the quarter-wave 
retarder plate relative to the beam splitter of the circular polarime¬ 
ter. The time to complete the cycle of eight subexposures for each 
observation accounted for 2h and 4 min. The achieved signal-to- 
noise ratio (5 /N) in the final Stokes I spectra summed over eight 
subexposures is rather low, accounting for a S IN = 14 during the 
night of 2012 July 15, and aS/N - 65 for the second night on 2012 
July 18. The spectra have a resolving power of about R = 115 000 
and cover the spectral range 3780-6910 A, with a small gap be¬ 
tween 5259 and 5337 A. The reduction and calibration of these 
spectra was performed using the HARPS data reduction software 
available at the FSO headquarter in Germany. 

The normalisation of the spectra to the continuum level con - 
sisted of several steps described in detail bv iHubrig et al.l l l2013bh . 
The Sto kes I and V parameters were derived followin g the ratio 
method jTinbergen & Rutteiil 1 19921 : IPonati et"^ Il997h . and null 
polarisation spectra were calculated by combining the subexpo¬ 
sures in such a way that polarisation cancels out. These steps e nsure 
that n o spurious signals are present in the obtained data (e.g. 

l2012h . 


Table 1. Magnetic field measurements of the Herbig Ae star PDS 2 us¬ 
ing the SVD method and the moment technique. All quoted errors are Icr 
uncertainties. 




SVD 

Moment technique 

HJD 

S/N 

S/Nsvd 

{Rz)sVD 

[G] 

<Bz>Fe 

IG] 

<Jlz)Fe,n 

IG] 

2456123.906 

74 

2420 

5±5 

23±25 

22±26 

2456126.904 

65 

2350 

33±5 

85±26 

18±28 


3 MAGNETIC EIELD MEASUREMENTS 

The software packages used by our group to st udy stell ar m agnetic 
fields are the moment technique developed bv lMath^ (e.g. ll99ih 
and the so-called multi-line Singular Value Decomposition (SVD) 
method for Stokes profile reconstruction recently introduced by 
ICarroll et alj ( |2012|) . Usually, before we apply the SVD method, 
we examine the presence of wavelength shifts between right- and 
left-hand side circularly polarized spectra (interpreted in terms of 
a longitudinal magnetic field (S^)) in a sample of clean unblended 
spectral lines using the moment technique. This technique allows us 
not only the estimation of the mean longitudinal magnetic field, but 
also to prove the presence of crossover effect and of the quadratic 
magnetic field. This information is highly important as, depending 
on the magnetic field geometry, even stars with rather weak longi¬ 
tudinal magnetic fields can exhibit strong crossover effects and kG 


quadratic fields (see e.g. 

Mathvslll995l:lLandstreet & Mathvsll2000l: 

iMathvs & Hubridll995l 

20061). In the studv of PDS 2, for each line 


in the selected sample of 44 unblended Fe ii lines, the measured 
shifts between the line profiles in the left- and right-hand circularly 
polarized HARPS spectra are used in a linear regression analysis 
in the AT versus d^geg d iagram, following the formalism discussed 
bv iMathvsl ( I1991 L[i 9^ . Our measurements <fiz) = 23 ± 25 G and 
(fiz) = 82 ± 26 G achieved on the first and second epochs, respec¬ 
tively, indicated the probable presence of a weak magnetic field on 
the second epoch. No crossover or mean quadratic magnetic field 
at a significance level of 3cr has been detected on either observing 
night, probably due to the rather low S/N of the observed spectra. 
The measured values for the mean longitudinal magnetic field 
are presented in the fifth column of Table[T] In the last column, we 
also present the results of the magnetic field measurements using 
the null polarisation spectra, labeled with n. Since no significant 
fields could be determined from the null spectra, we conclude that 
no noticeable spurious polarisation is present. 

The basic idea of the SVD method i s similar to the Principal 
Component Analysis (P CA) approach jMartmez Gonzalez et al.l 
120081 : ISemel et ^ l2009l) . where the similarity of the individual 
Stokes V profiles allows one to describe the most coherent and sys¬ 
tematic features present in all spectral line profiles as a projection 
onto a small number of eigenprofiles. 

As we mentioned above, the star PDS 2 was reported to show 
S Scuti-like pulsations on a time scale between 0.9 and 2.0 h. Since 
the full HARPS exposure times of about 2 h on both nights are of 
the same order, and pulsations are known to have an impact on the 
analysis of the presence of a magnetic field and its strength (e.g., 
ISchnerr et alj200^ : lHubrig et al.l2011bh . as a first step, we verified 
that no change in the line profile shape or radial velocity shifts are 
present in the obtained spectra. In Fig.[T]on the left side we present 
overplotted Stokes I profiles computed for the individual subex¬ 
posures recorded on both nights. On the right side of this figure we 
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Figure 1. Comparison of the SVD Stokes I profiles in the subexposures 
recorded on the nights of 2012 July 15 (top) and 18 (bottom). Left panels: 
Oveiplotted Stokes I profiles computed for the individual subexposures ob¬ 
tained with a time lapse of 15 min. Right panels: Differences between the 
Stokes I profiles computed for the individual subexposures and the average 
Stokes I profile. 
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Figure 2. /, V, and N SVD profiles obtained for PDS 2 on two different 
nights. The V and N profiles were expanded by a factor of 30 and shifted 
upwards for better visibility. The red (in the on-line version) dashed lines 
indicate the standard deviations for the V and N spectra. 


determined with SVD, and {Bft determined with the moment tech¬ 
nique from iron lines. 


show differences between the Stokes I profiles computed for the in¬ 
dividual subexposures and the average Stokes I profile. No impacf 
of pulsafions at a level higher than the spectral noise is detected in 
the Stokes / profiles. 

The resulting mean Stokes I, Stokes V, and Null profiles 
obtained by using the SVD method with 2094 metallic lines in 
the line mask are presented in Fig. |2] The line mask was con¬ 
structed using the VALD database (e.g. iKunka et al.l I 2 OO 0 I) and 
the respective stellar parameters of PDS 2. Due to the very low 
signal-to-noise at the blue and red ends of the HARPS spectra, the 
spectral lines used for the computation of the SVD profiles have 
been selected in the wavelength region from 4200 to 6200 A. The 
mean longitudinal magnetic field is estimated from the SVD re- 
con structed Stokes V and I using the center-of-gravity method (see 
e.g. lCarroll & Strassmeieil2014 ). A velocity range ±21 km s“* was 
adopted to determine the detection probability and the longitudinal 
magnetic field value. For both observations the null spectra appear 
flat, indicating the absence of spurious polarisation. No longitudi¬ 
nal magnetic field is detected on the first epoch on 2012 July 15, 
where we measure {Bft = 5 ± 5 G. A definite detection of a field, 
(Bz) = 33 ± 5 G, with a false alarm probability (FAP) smaller than 
10“®, is achieved on the second epoch three nights later. 

In Table [T] we list for the two observations the heliocentric 
Julian date, the 5 /N reached in the final Stokes / profiles, the S /N 
obtained in the SVD profile, the longitudinal magnetic field {Bft 


4 MAGNETIC FIELD MEASUREMENTS REPORTED IN 
PREVIOUS SPECTROPOLARIMETRIC STUDIES 

As the magnetic field strength and the magnetic field geometry 
in Herbig stars are poorly known - only about 20 late Herbig 
Be and Herbig Ae stars have been reported to possess large-scale 
organized magnetic fields - we compiled for them all magnetic 
field measurements reported in previous spectropolarimetric stud¬ 
ies. In Table |2] we present the {B^} values selected from previous 
low-resolution spectropolarimetric studies with FORS 1/2, where 
magnetic field measurements were carried out using the entire 
spectra, and high-resolution spectropolarimetric studies using the 
HARPS, ESPaDOnS, and Narval spectrographs. The rms longitu¬ 
dinal magnetic field, the rms standard error, and the reduced for 
th ese measurements have been computed following the equations 
of iBorra et al.Ul983l) : 
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Figure 3. Density distribution of the rms longitudinal magnetic field values 
(Col. 4 in Table |2) for the twenty late Herbig Be and Herbig Ae stars for 
which detections of a magnetic field were reported in the past. 


Tablel^lists for each star its name and spectral type, the num¬ 
ber of magnetic field measurements used to calculate the rms 
longitudinal magnetic field (Bz), the rms standard error, the reduced 
values, and finally the references to the measurements. 

In Fig. 0 we present the distribution of the rms longitudinal 
magnetic field strength for all late Herbig Be and Herbig Ae stars, 
for which detections of a magnetic field were reported in the past. 
The obtained density distribution of the rms longitudinal magnetic 
field values reveals that only very few stars have rms fields stronger 
than 200 G, and half of the sample possesses magnetic fields of 
about 100 G and less. 


5 DISCUSSION 

In this work, we present the first detection of a weak magnetic 
field in the Herbig Ae star PDS 2 using high-resolution spectropo- 
larimetric observations with HARPS. Previous observations using 
low-resoluti on spectropolarime try with FORS 1 indicated = 
103 ± 29 G jHubrig et al.l2009h . Apart from PDS 2, the authors an¬ 
nounced the detection of a weak longitudinal magnetic fields in six 
other Herbig Ae stars. Confirmation of the definite presence of the 
magnetic field based on accurate high-resolution spectropolarime¬ 
try with low uncertainties is still pending for these other six stars. 
Among these Herbig Ae/Be stars, three show longitudinal magnetic 
fields below 100 G, while for the remaining three stars the magnetic 
field ranges from 100 to 200 G. 

Magnetospheric accretion has been well established for 
T Tauri stars and depends on a strong ordered predominantly dipo¬ 
lar field channeling circumstellar disk material to the stellar sur¬ 
face via accretion streams. However, no convincing scenario on 
how magnetospheric accretion works in Her big Ae stars possessing 
rather weak magnetic fields currently exists. ICaulev & Johns-Krulll 
( I 2 OI 4 I) recently studied He i /H0830 morphology in a sample of 
56 Herbig Ae/Be stars. They suggest that Herbig Be stars do not ac¬ 
crete material from their inner disks in the same manner as T Tauri 
stars, while late Herbig Be and Herbig Ae stars show evidence for 
magnetospheric a ccretion. Further, due to their high rotation rates 
Im uzerolleir2004t) and weak magnetic fields, more compact magne¬ 
tospheres in Herbig stars are proposed. We note that since the rota¬ 


tion period of PDS 2 is currently unknown, it is not clear whether it 
might be a rapid rotator with a low inclination of the rotation axis. 

In fact, the density distribution of the rms longitudinal mag¬ 
netic field strength clearly indicates that the magnetic fields in Her¬ 
big Ae stars are by far weaker than those measured in their lower 
mass counterparts, the T Tauri stars, usually possessing kG mag¬ 
netic fields. As is shown in Fig. out of the sample of twenty 
magnetic late Herbig Be and Herbig Ae stars, only very few stars 
have rms fields stronger than 200 G, and half of the sample pos¬ 
sesses magnetic fields of about 100 G and less. The obtained results 
seem to support the existence of a different accretion mechanism 
mediating the mass flow from the disk to the intermediate-mass 
star, compared to that working in TTauri stars. Noteworthy, for the 
currently best studied Herbig Ae stars, HD 101412 and V380Ori 
with strong magnetic fields monitored over the rotation cycles, the 
obtained magnetic field models assuming a simple centered dipole 
indicate rather large obliquity of the magnetic axis to the rota¬ 
tion axis dHubrig et al.l201 Q : lAlecian et alj2009l) . The fact that the 
dipole axes are located close to the stellar equatorial plane is very 
intriguing in view of the generally assu med magnetospheric a ccre- 
tion scenario. As was shown in the past dRomanova et alj2003ll . the 
topology of the channeled accretion critically depends on the tilt 
angle between the rotation and the magnetic axis. For large inclina¬ 
tion angles //, many polar field lines would thread the inner region 
of the disk, while the closed lines cross the path of the disk matter, 
causing strong magnetic braking. Clearly, the qualitative picture is 
expected to be different if the magnetic field topology was proven 
as more complex than a simple centered dipole model. 

For a better understanding of the role of magnetic fields in 
Herbig Ae/Be stars, it is important to carry out additional highly 
accurate magnetic field measurements of a representative sample 
of Herbig Ae stars over their rotation periods using high quality, 
high-resolution polarimetric spectra. Using multi-epoch observa¬ 
tions, it will become possible to disclose the magnetic topology on 
the surface of Herbig Ae stars necessary to understand the interac¬ 
tion of the magnetic field with winds, accretion disks, convection, 
turbulence, and circulation. 
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Table 2. rms longitudinal magnetic field strength, rms standard error, and reduced values of late Herbig Be and Herbig Ae stars. N gives the number of 
measurements for the individual targets, separately for the low resolution spectrographs FORS 1 and FORS 2 and the high resolution spectrographs HARPS, 
ESPaDOnS, and Narval. 
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